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We employ the Boltzmann equation for describing hadron production from a quark-gluon plasma 
(QGP) in ultrarelativistic heavy-ion coUisions. We propose resonance formation in quark-antiquark 
scattering as the dominant meson-production channel, which, in particular, ensures that energy is 
conserved in the recombination process. This, in turn, facilitates a more controlled extension of ha- 
dronization to low transverse momenta (pt), and to address the experimentally observed transition 
from a hydrodynamic regime to constituent quark-number scaling (CQNS). Based on input dis- 
tributions for strange and charm quarks with azimuthal asymmetries, 112 (pt), characteristic for 
RHIC energies, we recover CQNS at sufficiently high pr, while at low pr a scaling with transverse 
kinetic energy is found, reminiscent to experiment. The dependence of the transition regime on 
microscopic QGP properties, i.e. resonance widths and Q-values in the q -I- g ^ Af process, is 
elucidated. 

I. INTRODUCTION 
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, Among the surprising experimental results in Au-Au collisions in the first years of operation of the Relativis- 

■ tic Heavy-Ion Collider (RHIC) were several unexpected features of hadron spectra in the intermediate transverse- 
momentum regime, pt ~ 2 — 6 GeV P, 0. Most notably, a large baryon-to-meson ratio of around one and the 

• so-called constituent quark-number scaling (CQNS) of the elliptic flow coefficient, V2{pt), have been found. These 
Cin' features are difficult to explain in hydrodynamic models 0, 0, Q which provide an excellent description of hadronic 
^ . observables in the \ow-pT region, pT < 2 GeV, but lack sufficient yield at intermediate pt- Quark coalescence models 
(]j ' (QCMs) have been used extensively and successfully to describe the baryon-over-meson enhancement and CQNS at 
f~i intermediate px [E 0) B B • The basic mechanism underlying QCMs is the recombination of constituent quarks 
into hadrons at the putative phase boundary (gluonic degrees of freedom are either neglected or converted "by hand" 
into q-q pairs). The constituent quarks are identified with the valence quarks inside hadrons, and, consequently, the 
hadron spectra directly reflect upon the partonic spectra and collective dynamics of the Quark-Gluon Plasma (QGP). 
T-H A major challenge is then to search for (and possibly construct) a unified description of hadron production at low and 
, intermediate px. In fact, as demonstrated in a recent PHENIX paper [ll|, a very general CQNS seems to emerge (in 
' centrality, system size and collision energy) when the scaling is applied to the kinetic energy of the produced hadrons. 
, One of the limitations of conventional QCMs is the collinear and instantaneous approximation of the hadron for- 
ty"} mation process which limits their applicability to intermediate pt where the inherent violation of energy conservation 
, is expected to be small. It is also not obvious how to recover an equilibrium limit in this framework which renders the 

■ connection to the lowp^ (hydrodynamic) regime less transparent. In this work, we propose an alternative realization 
of the recombination picture by formulating hadron production in terms of a Boltzmann transport equation. Our 
key dynamical ingredient is hadronic resonance formation via q-q "annihilation" . Their finite width not only enables 
to go beyond the collinear limit and explicitly conserve 4-momentum, but also introduces, in principle, a sensitivity 
to spectral properties of the QGP (e.g., widths and Q values in the fusion process). In the limit of large times, on 
which we will focus here, we recover an equilibrium formula which reproduces the thermal Boltzmann limit and can be 
compared to QCM results. We will constrain ourselves to the cases of charm and strange quarks hadronizing within 
a bulk medium where the input anti-/quark phase distributions are parametrized with collective properties inferred 
from RHIC data. We will compute both px spectra and V2{pt) for 0, D and J/ip mesons, and study their scahng 
properties in connection with the transition from low to intermediate px. Our approach bears some similarity with 
recent work in Ref. which, however, does not address the V2 problem and utilizes different interactions. 



II. BOLTZMANN EQUATION WITH QUARK-ANTIQUARK RESONANCES 

Our starting point is the Boltzmann equation for the meson phase-space distribution, /m, 

+ v- VfM){x,p) = -T/jp fM{x,p) + P{x,p) , (1) 

where p and x denote 3-momentum and position of the meson M, and F its width which in the present work 
is attributed to 2-body decays into quark and antiquark, M ^ q + q, and assumed to be constant. The factor 
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7p = EM{p)/m (m, Em{p)'- meson mass and energy), is due to Lorentz time dilation (see also Ref. The surface 

term on the left-hand-side of Eq. ([1]) is mostly relevant for high-pT^ particles, and vanishes in the infinite volume limit; 
since our objective here is not a quantitative description of data we neglect it in the following. The relation of the 
gain term, (3(x,p), to the underlying microscopic interaction becomes explicit upon integration over phase space, 

(fi p X ^, , [ (fi X (P pi P2 



^ ^ I !0^\-i f^iP^^) = / ^^^^6 fq{x,Pl)fq{x,P2) Cr(s) Vrel{pi,P2) (2) 



with (t(s) the cross section for the process q + q ^ M ai center-of-mass (CM) energy squared s = (p^'*'' + p'^-pY, 
p^^\ the 4-momenta of quark and antiquark, and /^^g their phase space distribution functions, normalized as Nqq = 

J ^J^i' fq,q{xTp)- Throughout this paper, quarks will be assumed to be zero-width quasi-particles with an effective 
mass Viiq which we treat as a parameter. The intrinsically classical nature of the Boltzmann equation warrants the 
use of classical distribution functions for all the particles. In addition, we will assume zero chemical potential for all 
quark species. For the cross section we employ a relativistic Breit-Wigner form, 

47r {Tmf 

^(«) -ff-p (g_^2)2 + (r^)2 ' (3) 

where g„ — gu I {9q9q) is a statistical weight given in terms of the spin (-color) degeneracy, gM igq,q)i of the meson 
(anti-/quark) ; k denotes the quark 3-momentum in the CM frame. With M ^ q + q being the only channel, it 
follows that Tin = Tout = T. Detailed balance requires the same F in the loss term on the right-hand-side of Eq. ([T]), 
thus ensuring a proper equilibrium limit with t = 1/F the pertinent relaxation time. This formulation conserves 
4-momentum and applies to all resonances M with masses close to or above the qq threshold, i.e., for positive Q value, 

Q = m — (niq + rriq) > 0. (4) 

If the 2^1 channel proceeds too far off-shell, i.e., Q < and F < \Q\ (e.g., for pions), other processes need to be 
considered, e.g., q + q M + g (which is possible in the present formalism by including the respective cross sections). 
Let us now elaborate the equilibrium limit of our approach, by imposing the stationarity condition, 

dNM I f (Px (Pp ^^eq 

dt ^'"^ J (271^3 



= -^l^-^ = - / -iTrrf'^fMhp + B. (5) 



Introducing the notation pM = Vi ^ P2iPrei — Pi ~ P2 into the gain term eq. ([2]), we find 

/d^ X d^p][j 
(27r)3 ^m{pm,x) ; /||(p, x) = g{pM, x)-/p/r (6) 

f d^Prel 

g{PM,x) = J -^j^^ fq{x,pM,Prel) fq{x,PM,Prel) Cr{s) Vrel{pM,Prel) ■ (7) 

Eq. ([6]) represents the large time limit of the Boltzmann equation and is the expression which comes closest to the 
conventional QCM formula. For hadronization times less or comparable to the relaxation time, r, the equilibrium 
limit will not be reached and a short time solution will be appropriate. In that case, the time variable enters explicitly 
into the final result, reflecting the dynamical nature of the Boltzmann equation. 

We have verified numerically that for F ^ Eq. ([6]) accurately recovers the standard Boltzmann distribution for 
a meson M at temperature T, if the constraint of a positive Q value is satisfied (for negative Q the 2^1 channel is 
inoperative). This shows that equilibration and energy conservation are closely related in our approach, constituting 
a significant improvement over previous QCMs. 



III. QGP FIREBALL, SPECTRAL PROPERTIES AND TRANSVERSE-MOMENTUM SPECTRA 

To compute meson spectra we have to specify the input q and q distributions (including their masses and collective 
properties), as well as the meson- resonance masses and widths. For an exploratory calculation, we will focus on a QGP 
fireball close to the expected hadronization temperature, Tc ~ 170 MeV (as in previous QCM studies). For simplicity, 
we will assume a constant (average) homogeneous cylindrical volume with collective transverse expansion velocity 
with linear radial profile, vt{tt) — Poj^ {Rt, Po'- surface radius and flow velocity), and use the quasi-equilibrium 
limit, Eq. ([6|). The quark distributions in the local (thermal) frame take the form 



fq{p, x) = exp [-jT{Eq -p- vt)/T] 



(8) 
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where 7t = (1 — v'^)~^/'^ and Eq = {rriq is the quark energy in the lab frame. In this section we evaluate 

Pt spectra for central ^/snn = 200 GeV Au-Au collisions, thus neglecting any azimuthal asymmetries. While the 
shape of the pr spectra is affected by f3o, the total multiplicities must be independent of /3o- Based on the fireball 
eigenvolume, V* — J (Px*, corresponding to a system at rest {(3q = 0), the volume element in the lab frame follows 
from the corresponding Lorentz-contraction factor (Pxiab = (Px* /^{x). We have verified that this leads to total 
quark multiplicities independent of /3o. The eigenvolume is fixed at V* ~ 1500fm'^ using a hadron resonance gas in 
chemical equilibrium at T ~ 170 MeV to reproduce the correct pion rapidity density dNT^/dy{y = 0) for 0-10% central 
collisions (including feed-down) [isj. Since thermal rapidity distributions for heavy particles are narrower than for 
thermal pions, the dN/dy{y — 0) values for the former have to be renormalized by a rapidity-width ratio, r*^/r^ 
(r^ = 1.8), to recover the proper equilibrium ratio for dNM /dy{y = 0). In this way, our calculated meson spectra are 
absolutely normalized. 

Based on Eq. the invariant meson pT spectrum takes the form 

Em^^ = "^^^^ = / ^ gip x)j /r (9) 

d^p d^pT dy (27r)3 J ^{x) ' ^ 

with g{p, x) defined in Eq. In the azimuthally symmetric case, 3 integrations can be performed analytically, leaving 
a 3-dimensional integral to be computed numerically. With a dependence on the transverse radius of the form (■^), 
the eigenvolume of the fireball factorizes and enters as an overall factor for pT spectra. A more complete description 
of Pt spectra should account for a hard component, which at RHIC energies is expected to become significant at 
PT > 2 3 GeV. Since we here focus on the conceptual aspects of our approach, we will defer the inclusion of a hard 
component to a future paper. 

Let us now specify the resonance parameters. The </> meson, which in the vacuum has a mass and width of 
rriff, ~ 1020 MeV and ~ 4 MeV, is expected to broaden substantially in hot and dense matter (cf., e.g., Ref. p^). 
while mass shifts are more controversial. There are indications from lattice QCD that the survives in the QGP at 
moderate temperatures above \\^ . We will assume for the its vacuum mass and a default width of = 50 MeV, 
in connection with a strange quark mass of rag = 400 MeV. The concept of D-meson resonances in the QGP has 
been implemented in Ref. [iq as a means to understand kinetic charm-quark equilibration at RHIC, with fair success 
in predicting [l7j the most recent data for (decay-) electron suppression and V2 [18]. Here we will assume _D-mesons 
with a mass and width of mo ~ 1.9 GeV and Yd — 100 MeV, and c- and u-quark masses of rric = 1.5 GeV and 
ruu = 350 MeV. For the J/i/' we employ mj/^, = 3.1 GeV, Fj/^ = 100 MeV. 

Fig. [T] illustrates our results for (left panel), D and J/ 4' (right panel) pT spectra compared to data in central 
y^'JNN = 200 GeV Au-Au. The surface flow value of (3o = 0.55 has been adjusted to best reproduce the (j) pT spectra; 
it is comparable to the value used in the conventional QCMs [7|. The agreement is fair up to the highest currently 
measured pt — 3 GeV. The spectra are rather insensitive to variations in F and Q (provided that F < Q), but this will 
be different for the elliptic flow discussed below. For J/ip and D spectra we have introduced a charm-quark fugacity, 
7c, which is necessary to match the total c-quark number to the one expected from hard production in primordial 
nucleon-nucleon collisions [1^, Note that Nj/^ oc 7^, with 7c subject to appreciable uncertainty in the charm 
cross section as well as due to in-medium effects on open-charm states 25], cf. also Ref. l26j . 



IV. ELLIPTIC FLOW AND SCALING PROPERTIES 



Azimuthal asymmetries in the particle's momentum distributions are defined via Fourier moments of order n, 

VniPr) {cos{rnp)){pT) , (10) 

where the average is over EdN/d'^p with d'^p/E = dydp'^dtp. At midrapidity, the odd moments vanish; in the following 
we will focus on the elliptic flow, V2{pt)- As discussed in the introduction, the experimental finding of an approximate 
CQNS relation for essentially all light hadrons at intermediate pt, 

V2{pt) = riQ vl{pT/nQ) (11) 

is naturally reproduced by QCMs and constitutes one of their successes (ng: number of constituent quarks in hadron 
h). As in QCMs, our objective is to study how a given input in the quark distributions manifests itself on the 
hadron level. We adopt a factorized ansatz for the quark u|, as used in previous coalescence studies 0,0, 

fq{p,x,^) = fq{p,x) ^{1 + 2 v^ipT) cos(2^)) (12) 

with fq{p,x) the distribution function of Eq. ([5]). Note that this "local" implementation of V2 (cell by cell in the 
fireball) neglects space-momentum correlations characteristic for hydrodynamic expansion (more realistic distributions 
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Pt (GeV) Pt (GeV) 

FIG. 1: Left panel: pr spectra for resonance formation at Tc — 170 MeV within our formalism, compared to data in central 
Au-Au at RHIC _19,;20',!.21]. Right panel: pr spectra for J/\& (lower curve) and D (upper curve) resonance formation, compared 
to PHENIX J/*" data in central Au-Au 22]; 'y^T = 170 MeV) = 7.8 has been taken from Ref. 



will be studied in forthcoming work). For the function V2{pt) we employ a simplistic ansatz, which, nevertheless, 
encodes the most importantphenomenological features. These are an essentially linear increase at low pT (as found 
in hydrodynamic models [1, 0, S HI] ) j ^^^^ a saturation at intermediate pr, requiring two parameters: the transition 
momentum, and the saturation value of U2, cf. also Refs. [2^.[29j. For light quarks (it, d, s) we choose p^* = 1.1 GeV 
and V2^^ = 7%, respectively, which gives a reasonable schematic representation of the experimentally found scaling 
function [ll[. Similar features are borne out of parton cascade models (3p| . [ Sll . For charm quarks, we take guidance 
from Langevin simulations for heavy quarks in a thermal background [13, [S^ • Also here, a plateau value of vf^*' = 7% 
is within the uncertainty of current data, while the transition momentum is shifted to a higher value, p^* ~ 2 GeV, 
which is compatible with the mass ordering effect in hydrodynamic calculations We have verified that our 

schematic parametrization is consistent with a parton-level KEt scaling for s and c quarks (characteristic for a 
hydrodynamically expanding parton phase), where KEt = — "t, denotes the transverse kinetic energy of a 
particle (m^ = \/ rn? + Pt)- "^2 saturation, implying deviation from hydrodynamics, signals, of course, incomplete 
thcrmalization and the transition to a kinetic regime. In Fig. [2] we summarize our results for the V2{pt) of (j) ^-nd 
J/^P mesons. The dependence of on F is very weak, changing by less than 5% when varying the (p width over 
the range F = 20 — 400 MeV. The dependence of V2^ on the Q value for the resonances is more pronounced as 
indicated by the different curves which have been obtained by varying the underlying quark masses (e.g., for the 0, 
Q = 220 MeV implies rus — (m^ — (5)/2=400 MeV). For any value of Q, CQNS, Eq. (fTT|) . is recovered in our approach 
at sufficiently high px where the production mechanism for an on-shell meson requires two constituent quarks with 
essentially collinear momenta, ~ ~ pt/2. Consequently, the meson will inherit the full azimuthal asymmetry 
imparted by its constituents. As is well known, for collinear production with V2 <C 1, using Eq. (|12p . the product fqfq 
recovers the scaling relation, Eq. (|lip . 

Deviations from coUinearity are expected to induce correction terms involving Q/px- Since 4-momentum is con- 
served in our approach these corrections can be quantified. In particular at low px there is no kinematical constraint 
enforcing coUinearity of the reaction. For Q — > 0, we accurately reproduce CQNS at all px, with the meson V2 
exhibiting twice the asymptotic value of u| reached at exactly twice the quark transition momentum. For increasing 
Q, however, the convergence to the limiting value is delayed to higher px, together with a reduction of the V2 at low 
px (resembling the effect of a larger mass). 

We finally return to the question of a universal CQNS. A recent analysis of the PHENIX collaboration [TJ] confirmed 
that CQNS for hadron px spectra is only satisfied at sufficiently high px, but the data exhibit a remarkably universal 
CQNS scaling when plotted vs. kinetic energy of the hadrons, KEx — rnx — m. At low px, KEx scaling is suggestive 
for hydrodynamic behavior [33| , followed by a transition to a kinetic regime where the V2 levels off. In Fig. [3] we 
summarize the CQNS properties as computed in our approach. The left panel displays the V2 for 0, D, and J/vI* 
when scaled in px- The 3 curves are different at low pr, converging at momenta as high as pr 5 GeV. However, 
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when the scahng is appUed in KEt (right panel of Fig. [3]), the curves essentiaUy coincide over the entire energy 
range, in quahtative agreement with experiment. While suggestive for the transition from the hydrodynamic to the 
kinetic regime, it remains to be scrutinized in how far this result depends on the uncertainties in the underlying 
parametrization of the quark W2's. As discussed above, the latter have been approximated by hydro-like behavior at 
low pt with a sharp transition to an asymptotic value of 7% at a momentum as estimated from transport calculations 
and data. An interesting point is that a KEt scaling implemented at the partonic level appears to manifest itself as 
KEt scaling at the meson level if the Q value in the q + q ^ M reaction is not too large (typically below 300 MeV, 
recall Fig. [2]). Another critical assumption is equality of the saturation value for wl for q = u, s,c, which, in turn, 
is suggestive for full c-quark collectivity in a strongly coupled QGP (sQGP). We finally note that, especially for the 
J/ip, the contribution from regeneration processes (as calculated here) is not expected to account for the yield at 
Pt ^ 4 GeV (cf. right panel of Fig. [1]), which implies deviations from CQNS at high momenta. 
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V. CONCLUSIONS 



We have proposed a reformulation of quark coalescence approaches based on a transport equation by implementing 
resonant quark-antiquark interactions to form mesons as the key microscopic mechanism. Our approach improves 
previous coalescence models by including explicit energy conservation in the meson production process and a well- 
defined thermal equilibrium limit. In the present exploratory study we have focused on a quasi-stationary scenario by 
considering the meson formation process in a mixed phase at constant temperature, T ~ T^. We recover the empirically 
found constituent quark-number scaling (CQNS) in the meson V2 at high momenta, while the improvements allow to 
extend the description to the low-momentum regime. We have quantified deviations from CQNS in px in terms of 
underlying spectral properties of the meson resonances, i.e., their widths and Q values. Larger values of the latter 
are found to shift the onset of CQNS to higher momenta. An interesting result of our study is that, for input (c- 
and S-) quark distributions satisfying kinetic-energy scaling, a KEt scaling is recovered at the meson level (0, £), 
J li^) provided the Q-values are positive and not too large. A universal KEt scaling of hadron V2 has recently been 
established experimentally [ll|. To consolidate our results, the sensitivity to the input V2{pt) at the quark level 
(including space-momentum correlations) needs to be scrutinized • Future work should also address the problem 
of baryon formation (e.g., within a quark-diquark picture [13])) an explicit time evolution (before and after T^) as 
well as hadron channels with negative Q values (e.g., for Goldstone bosons tt and K). While our framework does not 
address the hadronization problem of gluons, we believe that it could lead to useful insights into systematic features 
of hadron production over a rather wide range of momenta. 
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